Scaling behavior for a class of quantum phase transitions 
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We show that for quantum phase transitions with a single bosonic zero mode at the critical point, 
like the Dicke model and the Lipkin-Meshkov-Glick model, metric quantities such as fidelity, that 
is, the overlap between two ground states corresponding to two values Ai and A2 of the controlling 
parameter A, only depend on the ratio 77 = (Ai — A C )/(A2 — Ac), where A = A c at the critical 
point. Such scaling property is valid also for time-dependent quantities such as the Loschmidt echo, 
provided time is measured in units of the inverse frequency of the critical mode. 
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I. INTRODUCTION 

Due to the quantum nature, a Quantum Phase Tran- 
sition (QPT) defined as a drastic change of funda- 
mental properties of ground states, may have properties 
quite different from thermal phase transitions. This point 
has received much attention in recent years. In partic- 
ular, some concepts and quantities in the field of quan- 
tum information, e.g., entanglement and fidelity, have 
been found quite useful in characterizing the occurrence 
of a QPT. For example, the overlap between two ground 
states corresponding to two nearby values Ai and A2 of 
the controlling parameter A, 



i P (Ai,A 2 ) = |(0 Al |0 A2 )|, 



(1) 



has been proposed as a probe of quantum criticality 
the dramatic change of the wave function at a QPT im- 
plies a decrease of the overlap L p in the neighborhood of 
the critical point. Hence, the fidelity L v can be used to 
detect the occurrence of a QPT (see [2h1C| and references 
therein). 

Similarly to thermal phase transitions, an important 
aspect of QPTs is the dependence of relevant quantities 
on the controlling parameter A. For example, the charac- 
teristic energy scale usually takes the form |A — A c | v , with 
A c indicating the critical point and (p > a critical expo- 
nent. This is for the case in which only one value of the 
parameter A is of relevance. For quantities like fidelity, 
which instead depend on two values Ai and A2 of the con- 
trolling parameter, one should understand whether their 
behavior in the critical region encodes universal proper- 
ties about a QPT. With regard to the fidelity, the ques- 
tion is whether its drop near a QPT can be used not only 
to detect the QPT itself but also to determine the criti- 
cal exponents [l0| . In this context, understanding scaling 
properties is relevant. 



'Electronic address: wgwang@ustc.edu.cn 



When studying fidelity, either as the overlap (JXJ) of 
ground states or as the survival probability of an initial 
state prepared in the ground state |0a 2 ) of Hamiltonian 
H (A2) and evolved under a different Hamiltonian H(Xi), 
two values Ai and A2 of the controlling parameter are 
involved. For the sake of clearness, in what follows, we 
use the name fidelity for the former quantity, L p (Xi, A2), 
and use Loschmidt Echo (LE) for the latter one [ll( • An 
interesting question is about the dependence of fidelity 
and LE on Ai and A2. In the case that fidelity L p (Ai, A2) 
goes to zero when Ai approaches A c for a fixed A2, it is 
clear that for each A 2 there exists a value A' : such that 
L p {\\,\2) = L p (A' 1 ,A' 2 ). This implies that the relative 
positions of Ai and A2 with respect to A c , rather than 
their exact positions, play the crucial role. Hence, the 
question arises of whether fidelity may be invariant un- 
der rescaling of the controlling parameter. 

In this paper, we show that for QPTs possessing 
only one bosonic zero mode at the critical point, met- 
ric quantities like fidelity only depend on the ratio i] = 
(Ai — A C )/(A2 —A c ). That is, these physical quantities are 
invariant under linear rescaling of the controlling param- 
eter with respect to the critical point. We also show that 
such scaling property is valid for time-dependent quanti- 
ties such as the LE, provided time is measured in units 
of the inverse frequency of the critical mode. The class of 
QPTs possessing such features includes important phys- 
ical models, like the Dicke [llj] and the Lipkin-Meshkov- 
Glick (LMG) model H3. 

The article is organized as follows. In Sec. [TT| we dis- 
cuss our scaling argument for static metric quantities 
like fidelity, extending in Sec. IIIII such scaling to time- 
depending quantities like the LE. The scaling for time- 
dependent quantities is then illustrated by means of the 
semicassical theory. We then illustrate the fidelity and 
LE scaling in two relevant physical models, the Dicke 
model (Sec.|TV| and the LMG model (Sec.©. We finish 
with concluding remarks in Sec. I VI I 
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II. SCALING ARGUMENT 

As only the lowest energy levels are concerned close 
to the critical point, wc assume that the Hamiltonian 
describing a QPT can be approximately written in terms 
of n harmonic oscillators: 



^(A) = ^e i (A)cJ(A)c i (A), 



(2) 



where c\(X) and Ci(X) are bosonic creation and annihila- 
tion operators for the i-th mode. The ground state |0 A ) 
for the parameter A is defined by Ci(A)|0 A ) = 0. For two 
parameter values Ai and A2, one may write 



n 

5t(A 1 ) = ^[^4(A 2 ) + Q ij c i (A 2 



(3) 



where Py and Qij are functions of Ai and A2 , with Pjj = 
5ij and Qy = for Ai = A 2 . Wc discuss the case in 
which there is only one zero mode at the critical point: 
ei(A) ~ I A — A c | v for A close to A c , while e,(A c ) 7^ for 
i In this case, Eq. © reduces to: 



cI(Ai) = Pnc t 1 (A 2 ) + Q 11 ci(A 2 ), 



(4) 



with a corresponding expression for ci(Ai). From the 
bosonic commutation relations it follows that |Pn| 2 — 
|Qn| 2 = 1- Let us write explicitly the phases of Pn 
and Q n , as P n = |Pn|e^+H Qn = \Q xl \e^~ e ^ . 
In the representation of P(A 2 ), the change of the pair 
(ci(A 2 ),c|(A 2 )) to (e" ie ''ci(A 2 ),e ie ''c t 1 (A 2 )) docs not bring 
any change to the physics, hence, the phase 9 r can be ab- 
sorbed by ci(A 2 ) and c{(A 2 ). The phase 9 C is the relative 
phase between the pair of operators (ci(A 2 ), c{(A 2 )) at 
A 2 and the pair (ci(Ai), c\(Xi)) at Ai, which generates 
relative phases between the set of basis states |« Al ) = 
( C t 1 (A 1 ))™|0 Al )|/V^! and \mx 2 ) = (c\(X 2 )) m \0 X2 )\/Vri.. 

Let us consider a physical quantity A depending on 
two values Ai and A 2 of the controlling parameter (for 
instance, A might be the fidelity), written in the vicin- 
ity of the critical point as a function of the annihilation 
operators for the zero mode: 

A = (0 Ai |i(c 1 (A 1 ),c 1 (A 2 ))|0 Aj ), (*,j = 1,2). (5) 

In what follows, we focus on quantities A that do not 
depend on the phase 9 C [Til ]. In particular, metric quan- 
tities, like Cnm = \{n\Am\ 2 )\, belong to this class. Such 

00 



quantities include, for instance, the fidelity L T 



and the participation ratio x °f an eigenstate of e.g. 
P(A 2 ), \mx 2 ) = E, l ("Ail"iA 2 )|n-A 1 ), with respect to the 
basis of the eigenstatcs of H(Xi); by definition, \ = 

!/E„ K^aJ^)! 4 = l/En C L' In thc stud Y of these 
quantities, we can take 6 C = 0. Then, since the phase 
8 r can be absorbed by ci(A 2 ) and c|(A 2 ), Pn and Qn 
are just their absolute values, with Pn > 1. Using the 



ground state definition ci(Ai)|0 Al ) = 0, Eq. (|4]), and the 
expansion |0 Al ) = J2 m ( mx 2 |0ai)|w A2 ), we can express 
|0 Al ) as a function of Pn, Qn, c|(A 2 ), and |0 Aa ). After 
inserting the obtained expression for |0 Al ) into ([5]), wc 
find that A is a function of Pn and Qn only. 

The dependence of Pn on Ai and A 2 can be written 
as Pn = P(AAi, AA 2 ), where AA, = A 4 - A c (i = 1,2). 
We study Ai and A 2 belonging to the same phase, so 
that t) = AAi/AA 2 > 0. We assume, as it is natu- 
ral for a QPT with an infinitely-degenerate zero mode 
at the critical point, that P(AAi,AA 2 ) goes to infinity 
in the limit AAi — > with AA 2 7^ 0, as well as in the 
limit AA 2 — > with AAi 7^ 0. For the sake of sim- 
plicity, we assume that F is a monotonic function of 
AA 2 , when AA 2 changes from a given AAi to fl5j |. 
Then, given P(AAi,AA 2 ) = d, for any X[ there must 
exist a A 2 such that P(AA' 1 ,AA' 2 ) = d. This implies 
that there exists a function AA 2 = g(AAi,d), such that 
P(AAi, g(AAi, d)) = d for any AAi, hence, 



8F/dAX 1 + {dF/dAX 2 )(dg/dAX 1 ) = 0. 



(6) 



For a given d and a sufficiently small AAi, Tay- 
lor expansion reads g = go + g'AXi, where g' — 
dg/dAXi with d fixed. We recall that in the limit 
AAi — > 0, P(AAi, AA 2 ) goes to infinity if AA 2 is non- 
zero. Hence, for any given d, AA 2 must go to zero 
in order that limA Al ^o P(AAi, AA 2 ) = d. That is, 
lim.A Al ^o <?(AAi, d) = 0, as a result, go — 0. There- 
fore, when the above Taylor expansion works for all fixed 
values of d, we have dg/dAXi = g/AX\ = AA 2 /AAi. 
Substituting this result into Eq. ((5]), we find 



AAi OF/dAXx + AA 2 8F/dAX 2 = 0. 



(7) 



This equation has the solution F = P(ln AAi — In AA 2 ) = 
P(ln?/). Hence, Pn is a function of 77. Then, due to 
the relation P n — Q 2 ^ = 1, Qn is also a function of 
77. Since the quantity A is a function of Pn and Qn, 
we can conclude that A only depends on the ratio 77 = 
(Ai - A C )/(A 2 - A c ) 



III. TIME-DEPENDENT QUANTITIES 

We now consider time-dependent metric quantities 
A(t), with the dynamics described by Hamiltonian ©. 
Since the frequencies 0Ji(X) = e,(A)/?i depend on A, A(t) 
usually cannot be a function of 77 only. However, for sys- 
tems with a single zero mode at the critical point, the 
77-scaling still applies, provided time is rescaled: t — > r = 

LJl(X)t. 

As an illustration we show in the following that in the 
vicinity of a QPT with a single zero mode at the critical 
point the decay of the quantum Loschmidt echo depends 
only on the scaling parameter 77 and the rescaled time r. 
The LE gives a measure for thc stability of the quantum 
motion under slight variation of the Hamiltonian (l6l - [l8| . 
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It is defined by M^t) = \m(t)\ 2 , where 

m(t) = (^ a \cxp(iH{X 2 )t/h)cxp(-iH{X 1 )t/h)\^ ). (8) 

Here, H{Xi) = H(X 2 ) + eV, with e = A x - A 2 . Exten- 
sive investigations have been performed in recent years, 
to understand the decaying behavior of the LE in dif- 
ferent regimes, depending on the chaotic or integrable 
nature of the dynamics, on the system's dimensional- 
ity, and on the perturbation strength (see Refs. [l9l - [3^ 
and references therein). Furthermore recent investiga- 
tions have shown that the LE may be employed to char- 
acterize QPTs, since it exhibits extra-fast decay in the 
vicinity of critical points [H|-[39| . 

Here we consider a system initially prepared in the 
ground state |0a 2 ) of H(X 2 ). Then the LE is in fact the 
survival probability, 



M L (t) = |(0 A2 | e -^ Al )*/ ri |0 A2 )| 2 



(9) 



In the critical region, H(X 2 ) represents a harmonic oscil- 
lator and, therefore, its ground state can be written as a 
Gaussian wave packet. As shown in Ref. [3l|, when the 
classical motion is periodic with a period T p , scmiclassi- 
cal theory predicts that for t > T p the LE has an initial 
Gaussian decay followed by a power law decay. Indeed, 
to a second-order term of perturbation expansion, 

M L (t) ~ 6 (1 + £h 2 )- 1/2 e- rt2/{1+et2 \ (10) 



where bn 



«o ~ i, r = (^) 2 A e 

the derivatives evaluated at the center p$ of the initial 
Gaussian wave packet. Here, U = J^ F V dt and W is 
a measure of the width of the initial Gaussian packet in 
the momentum space. It is seen that Ml has a Gaussian 
decay e~ rt " for short times and a l/£t decay for long 
times. 

Let us consider the case in which H has, close to the 
critical point, two lowest energy relevant modes, of which 
only the first one has zero frequency at the critical point . 
The first mode moves very slowly, so that effectively the 
system moves periodically with the frequency of the sec- 
ond mode. That is, for times much shorter than the 
period T± of the first mode (which diverges at the criti- 
cal point), the classical motion is approximately periodic 
with the period of the second mode. Hence, the period 
T p = 2w / u 2 (Xi) . Note that for long times the first mode 
dominates and the LE oscillates with a period related to 
T\ = 2ir/u>i(\i) (detailed later). Therefore, the above 
semiclassical prediction works within times longer than 
T p and shorter than T\. 

Starting from the classical expression of the Hamilto- 



1 eW 
I 2ft 



d 2 U\ 



with 



we obtain 



i & 2 / H(X 1 )-H(X 2 ) ^ 

T~ -r^r 7—^ r 



wi(Ai)' dpi ei(Ai 
&_/Ix{Xl) ei (A 2 )Jx(A 2 ) , AH 
dp 2 



ei(Ai) 



ei(Ai) 



(12) 
(13) 



where AH 2 = ^2 (^1)^2(^1) — -^2 ^2)^2^2)- The second 
mode has no singularity at A c , hence from Taylor ex- 
pansion of AH2 we obtain AH2 ~ (Ai — A2). Hence, 
^fjy ~ |1 - 77— 1 ]'^|Ai - X 2 \ 1 ~ v - Therefore, when X 1 is 

sufficiently close to A2, e ^ffi\ can be neglected for tp < 1 
and 77 7^ 0. Then, since ii(Ai) and Zi(A2) have no singu- 
larity at A c and — we find that in the very 

neighborhood of A c , £t ~ F(rj)r. Similarly, Ft 2 can be 
written as G(rf)T 2 . We can therefore conclude that 



M L (t) ~ 6o(l + (F(r ? )r) 2 )- 1 / 2 e- G (") T2 /(i+(F(r ; )r) 2 ) 



is a function of r\ and of the rescaled time r. 



(14) 



IV. SCALING FOR THE DICKE MODEL 

This model [l2j provides a physically significant ex- 
ample of our scaling behavior. It describes the interac- 
tion between a single bosonic mode and a collection of N 
two-level atoms and finds applications in quantum optics, 
condensed matter physics and quantum information. In 
terms of the collective operator J for the N atoms, the 
Dicke Hamiltonian is written as (hereafter we take H = 1) 

H(X) = w J z +u;a t a + (A/\/AO(a t + a)(i+ + J_). (15) 

In the thermodynamic limit N — > 00, the system un- 
dergoes a QPT at A c = ^^/HJuJo, with a normal phase 
for A < A c and a super-radiant phase for A > A r . The 
Hamiltonian can be diagonalized in this limit [l2[, tak- 
ing, up to a constant energy term, the form ([2]), with 
n = 2 modes. In the normal phase, the energies of the 
two harmonic oscillators read 



1/2 



ei.a(A) = { \ [(u 2 +u 2 )± yVo - w 2 ) 2 + 16A 2 ^ ] } 

(16) 

ordered so that ei(A) < e 2 (A). It is seen that ei(A) = 
for A = A c , hence, the ground level of H(X C ) is infinitely 
degenerate and the system undergoes a QPT at A c . On 
the other hand, 62(A) 7^ at the critical point. In the 
super-radiant phase, the Dicke Hamiltonian can still be 
diagonalized in the thermodynamic limit, resulting in a 
two-mode form, with the energies 



H(X) = wi(A)/i(A) + w 2 (A)/ 2 (A). 



(11) 



e li2 (A) = {\[. 2 + f 2 ± ^{f 2 - ^-) 2 + 4^]} 1/2 , 

(17) 

where /1 = loluq/AX 2 and ei(A) < e 2 (A). It is easy to see 
that ei(A) = and e 2 (A) ^ for A = A c . Thus, the 
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FIG. 1: (Color online) Dependence of the fidelity L p on the 
scaling parameter 77 for different values of Ai and A2. Note 
that the values of 77 on both sides of the zero point are positive. 
Symbols on the left hand side of 77 = represent fidelity in 
the normal phase, with the superscript N of A2 standing for 
normal phase (open symbols). Symbols on the right hand 
side of 77 = are for the super-radiant phase (superscript S of 
A2, full symbols). Data are in agreement with the analytical 
result of Eq. ()18|) (solid curve). 



ground level of H{\ c ) is also infinitely degenerate (and 
with a single zero mode) at the critical point from the 
super-radiant-phase side. 

As can be shown both analytically and numerically the 
fidelity in this model uniquely depends on the scaling 
parameter 77. Indeed, in both phases, it has been found 
(see Aappendix [Af that 



Lp — 



V2, 



>rft 



(18) 



with the critical exponent ip = 1/2 (o>i(A) ~ |A — 
A c | 1/2 ) S3]. For 77 < 1, L p oc rf /4 . The analytical result 
(|18[) is in agreement with numerical simulations shown in 
Fig. [T] data for different values of Ai and A2 collapse on 
a single universal curve. 

In the derivation of Eq. (|T8|) . we have used Eq. (IT6l) . 
which is obtained by diagonalizing an effective form of 
the exact Dickc Hamiltonian (|15p . given by Eq. @, with 
71 = 2 modes. It is therefore important to assess the valid- 
ity of the effective Hamiltonian in computing the fidelity 
L p . The effective Hamiltonian leads to an error inversely 
proportional to N, where N is the number of atoms. For 
those quantities, for which there exist some contributions 
proportional to N, the effective Hamiltonian gives poor 
predictions j4ll - l43l |. However, for quantities like the fi- 
delity L p of ground states, there is no such contribution. 
Hence, the effective Hamiltonian is expected to correctly 
describe the behavior of L p in the critical region. To sub- 
stantiate such expectation, we have compared the pre- 
diction L p (Ai,A2) of Eq. flTSl) with Lp(\ 1} \2), which is 
the corresponding fidelity numerically computed by di- 



FIG. 2: The quantity lnD versus IniV in the Dicke model, 
with parameters ojq = u = 1, Ai = 0.495, and A2 = 0.45 
(A c = 0.5). 



rect diagonalization of the exact Hamiltonian (1151) in a 
truncated Hilbert space. The truncated Hilbert space is 
obtained for a finite number N of the atoms and by tak- 
ing the lowest N levels of the bosonic mode. We have 
studied the variation of D with N, where 



D 



|i^(Ai,A 2 )- J L p (Ai,A 2 )|. 



(19) 



With increasing N, the quantity D exhibits a decay faster 
than power law (see Fig. [2]) and slower than exponential. 
Therefore, it is reasonable to expect that the effective 
Dicke Hamiltonian provides the correct physical picture 
when computing the fidelity of ground states in the large 
N limit. 

Next, we discuss the LE. For the Dicke model, it can be 
analytically proved that the LE is an oscillating function 
of time with period T = tz/u)\(\\) = Ti/2. This period 
diverges when Ai approaches the critical point A c and for 
times shorter than T/2 the LE decays according to the 
above semiclassical prediction. Indeed, numerical simu- 
lations in Fig. [3] show that the LE is a function of 77 and 
of the rescaled time r = ivi(Xi)t. Moreover, in the super- 
radiant phase the LE decays in the same manner as in 
the normal phase. Finally, we have studied the minimum 
value of the LE, denoted by M p , as a function of Ai and 
A 2 . Since this quantity is time- independent, according 
to previous scaling arguments it should be a function of 
the ratio 77 only. Such expectation is confirmed by our 
numerical simulations (see the inset of Fig. 



V. SCALING FOR THE LMG MODEL 

In the two-orbital Lipkin-Mcshkov-Glick (LMG) model 
for N interacting particles, in terms of the total spin 
operator for its collective motion, S a (a = x,y,z), the 
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FIG. 3: (Color online) Dependence of the Loschmidt echo Ml 
on the rescaled time r = u>i(Ai)£, for various values of Ai and 
A2, with different symbols representing different pairs (Ai, A2). 
The three curves correspond, from top to bottom, to r\ = 
1(T 2 , 1CT 3 , and 10" 4 . An initial Gaussian decay is followed 
by 1/t decay, as predicted in Eq. (|f4[) . Inset: Dependence of 
the minimum value M p of Ml on r\ for various values of Ai 
and A2 (open symbols stand for the normal phase and solid 
ones for the super-radiant phase) . The fitting curves are given 
by M p = 2y/rj/(l + rj) in both phases. 



Hamiltonian can be written as, 

(7, h) = ~(S 2 X + 7S 2 ) - 2hS z + (1 + 7 )/2. (20) 

As shown in Ref. [H|, in the thermodynamic limit, mak- 
ing use of the Holstein-Primakoff transformation and of 
a standard Bogoliubov transformation, the Hamiltonian 
can be diagonalized, 



H(-y,h) = Aa^a 



(21) 



where 



1-7 



A = 2\(h- l)(ft- 7 )] 1/2 , tanh9 = — '■ — (22) 

2h — 1 — 7 



for h > 1, 



A = 2[(l-/ l 2 )(l- 7 )] 1 / 2 , tanh6 = 



h 2 



7 



(23) 



for h < 1, and flg and ae are bosonic creation and anni- 
hilation operators. 

Equations (|2"2"T) and (f2"3")l show that when h approaches 
1 from both sides, A — s- 0. This implies that the sys- 
tem undergoes a quantum phase transition at the critical 
point h c = 1. The phase with ft, > 1 is usually called the 
symmetric phase and the phase with h < 1 the broken 
phase [H[. 

In the thermodynamic limit, the same scaling law as 
in Eq. (|18p can be derived analytically in the vicinity of 
the critical point. In fact, for a fixed 7 the ground state 



|0)e 2 for h = l%2 has the following expansion on the basis 
|n)ei of the eigenstates of Hamiltonian (f2"T|) at h = hi 



\2n- 1)!! 



tanh™ 



83-61 



(2n)!! 

where C is a normalization constant 

2 Ae 2 -9x 



c 



1 - tanlr 



-1/2 



|2n) ei (24) 



(25) 



Then, it is ready to find that 



L p (hi,h 2 ) 



1 — tanh 



2/62-6! 



1/4 



(26) 



In the vicinity of the critical point h c = 1, the right hand 
side of Eq. (f2l)f can be simplified further. In fact, in the 
symmetric phase, from Eq. (|22|) . one obtains, up to terms 
of higher order in h — h c , 



This gives 



tanh 



tanh 



9 



e 2 -6i 



1 



1-7 



1/2 



r? 1 / 2 - 1 
v i/2 +1 > 



(27) 



(28) 



where 77 = (/ii — l)/(/i2 — !)• After inserting ([28]) into 
(|26|) . we obtain the same expression (|18|) for L p as for the 
Dicke model. 

For the LE, making use of analytical results in the 
symmetric phase [45[ and its generalization to the broken 
phase, similar scaling behaviors as shown in Fig. [3] have 
also been found. 



VI. CONCLUSIONS 

To summarize, we have proved a scaling property for 
time-independent metric quantities such as the fidelity 
and the participation ratio. The scaling is valid for mod- 
els like the Dicke model and the LMG model, whose 
QPT can be described in terms of a single bosonic zero 
mode. Moreover, also time-dependent quantities such as 
the Loschmidt echo, exhibit the same scaling provided 
time is measured in units of the inverse frequency of the 
critical mode. 

Our scaling arguments showing ^-dependence of static 
quantities can be generalized to the cases of more than 
one zero mode, provided appropriate new restrictions in 
the coefficients Pij and Qij in © are introduced. On 
the other hand, the scaling for time-dependent quantities 
cannot be extended in a straightforward way to the case 
of more than one zero mode, when the corresponding 
frequencies tJi(X) have different scaling behaviors and, in 
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contrast to the case of a single zero mode, there is no 
natural rescaling of time by means of a single frequency. 

Finally, we remark that in our theory we compute the 
fidelity in the thermodynamic limit. The obtained results 
imply that many-body systems whose QPT can be de- 
scribed in terms of a single bosonic zero mode do not ex- 
hibit the Anderson orthogonality catastrophe [|||. That 
is to say, the ground states corresponding to two nearby 
values of the controlling parameter are not orthogonal 
at the thermodynamic limit, provided these two values 
belong to the same phase. 
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7 = — arctan[4A v / woJo/( w2 + w o)]- 
It is straightforward to verify the following relations: 

detA x = e$e%, (A2) 
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and 



WW, PQ, and QW acknowledge support by the 
Natural Science Foundation of China under Grant 
Nos. 10775123 and 10975123, the National Fundamental 
Research Programme of China Grant No.2007CB925200, 
and 'Boshidian' Foundation of the Ministry of Education 
of China. GB and GC acknowledge support by MIUR- 
PRIN 2008 and by Regione Lombardia. 



~ A 2 A 2 
O-i Co 



--2 
2 



(A4) 



In the normal phase of the Dicke model, energies e\ 2 
are given by Eq. ([To]) . In the neighborhood of the critical 
point A c , from Eq. (|16[) we get, up to terms of higher 
order in \ r — A, 



Appendix A: Derivation of Eq. (|18|) 

In this section, we give a derivation of Eq. (|18|) in the 
main text for the fidelity of two ground states at Ai and 
A2. For this purpose, we use the following expression of 
the fidelity given in Ref. 0], 



L, 



2{[dctA A2 ]/[dctA Al ]} 1/4 
[det(l + A^A X2 )]V2 ' 



(Al) 



where A x = U~ 1 MxU, M x = diag[e£,e£], and U is an 
orthogonal matrix, 



U 



c —s 
s c 



8A C (A C - A)ww 



1/2 



(A5) 



Then, we have e^/e^ 1 = (l/?y) 1/2 . Using again Eq. (|A3|) . 
we obtain 



det Ay 



det A\ 1 
Tt{AllAx 2 ) = y/ljv + l, 




(A6) 
(A7) 



where e^/e^ 1 = 1 has been used in the vicinity of 
the critical point. Substituting the above results into 
Eq. (|A1[) . one finds Eq. (fT8l . By the same method, the 
same expression of the fidelity can be obtained in the 
super-radiant phase. 
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